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The small, about 30% magnetoresistance at the onset of full spin polarization in the 2D electron 
system in a modulation-doped Si/SiGe quantum well gives evidence that it is the remote doping that 
determines the transport scattering time. Measurements of the mobility in this strongly-interacting 
electron system with remote-doping scattering allow us to arrive at a conclusion that the Hubbard 
form underestimates the local field corrections by about a factor of 2. 
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Much interest has been attracted recently by the be- 
havior of two-dimensional (2D) electron systems in a par- 
allel magnetic field. The resistance of a 2D electron sys- 
tem in Si MOSFETs was found to rise with parallel field 
saturating to a constant value above a critical magnetic 
field Be d] which corresponds to the onset of full spin po- 
larization of the electron system |2l . The spin origin of 
the effect is consistent with its insensitivity to the direc- 
tion of the in-plane field P|. At low electron densities, the 
experimental resistance ratio is equal to R(Bc) / R{Q) ~ 4 
as long as the system remains metallic [2j, which is in 
agreement with the calculation Q . The effect was used to 
study the spin susceptibility of the least-disordered elec- 
tron system as well as the local moments in the band- 
tail in more-disordered Si MOSFETs |i,l3|- In contrast, 
the 2D carrier system in GaAs/AlGaAs heterostructures 
is relatively thick so that the orbital effects become im- 
portant and give rise to an enhancement of the effec- 
tive mass in parallel magnetic fields As a result, 
in GaAs there are two noteworthy distinctions (i) 
above Be, the resistance keeps on increasing less steeply 
with no sign of saturation; and (ii) the magnetoresistance 
is strongly anisotropic depending upon the relative ori- 
entation of the in- plane magnetic field and the current. 

The 2D carrier systems in Si MOSFETs and high- 
mobility GaAs/AlGaAs heterostructures are similar in 
that at low carrier densities, the transport scattering time 
is determined by charged impurities near the 2D system. 
In the former system charged impurities are located at 
the Si/SiOo interface resulting in dominant large-angle 
scattering whereas in the latter they are homoge- 
neous background doping for both 2D electrons 12] and 
2D holes 13] . The case of a 2D carrier system with a 
finite spacer that is remarkable by remote-doping scat- 
tering, i.e., dominant small-angle scattering, is opposite. 
Although it may seem simple to realize such a 2D sys- 
tem, that kind of scattering has not been unequivocally 
established in any 2D system studied so far. Recently, 



it has been predicted that for remote doping, the above 
resistance ratio should be equal to R{Bc)/R{0) ~ 1.2 

In this paper, we report measurements of the resis- 
tance of the 2D electron system in a modulation-doped 
Si/SiGe quantum well in parallel magnetic fields. Being 
very similar to (lOO)-Si MOSFETs, this electron system 
is different by the presence of a spacer and the 1.5 times 
larger dielectric constant so that the same strength of 
electron-electron interactions in Si/SiGe can be expected 
at 2.2 times lower electron densities. The fact that the 
observed magnetoresistance is small shows that the trans- 
port scattering time in our samples is determined by re- 
mote doping. In this regime, the magnetoresistance is 
sensitive to a small amount of the charged residual im- 
purities at the Si/SiGe interface (< 1% of the electron 
density), which varies in different runs depending on ex- 
ternal perturbations such as cooling and illumination of 
the sample. Using the dependence of the mobility on 
electron density in this strongly-interacting system with 
remote-doping scattering, we extract the form of the local 
field corrections (LFC). From comparison of the exper- 
imental data and the model calculations it follows that 
account should be taken of both kinematic and dynamic 
correlations. 

Samples were grown by molecular beam epitaxy on 
(OOl)-Si substrates. A 15 nm thick Si channel was de- 
posited on a strain-relaxed buffer layer which consists of 
a 600 nm Sio.yGeo.s buffer on the top of a 2.5 /zm thick 
Sii_a;Gea; layer with compositional grading. The channel 
was capped by 14 nm Sio.rGeo.a spacer, 12 nm P-doped 
layer, and 27.5 nm cap layer, covered by 10 nm of Si. 
The samples were arranged in a standard Hall bar geom- 
etry. The channel of the transistor had a uniform width of 
20 /im between the ohmic source and drain contacts. The 
four voltage probes and the source/drain electrode were 
formed by local implantation of phosphorous ions. Sub- 
sequent activation and recrystallization was performed 
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FIG. 1: Zero-field mobiUty as a function of electron density 
at a temperature of 30 mK for two samples (open symbols) 
along with the data of Ref. [l^ (filled squares). Also shown 
is the result of the calculation with (solid lines) and without 
LFC (dashed line) neglecting the numerical factor. 



at 560° C for 30 min in an N2 ambient. Afterwards the 
contact pads were metallized with 25 nm Ti and 100 nm 
Au. In order to reduce the contact resistance, the sample 
was annealed at 400° C for 60 s in an N2 atmosphere. Fi- 
nally, the Schottky gate electrode comprising 30/100 nm 
Ni/Au was defined by a hft off process. 

The sample was placed in the mixing chamber of a 
dilution refrigerator with a base temperature of 30 mK. 
To create the mobile carriers in the 2D system, the sam- 
ple was illuminated with a light-emitting diode until the 
resistance saturated. After the diode was switched off, 
the sample state did not change in the run. The resis- 
tance was measured using a standard four-terminal lock- 
in technique at a frequency of 15 Hz in magnetic fields up 
to 14 T. Excitation current through the device was kept 
low enough (< 10 nA) to ensure that measurements were 
taken in the linear regime of response. To change the 
sample position in the mixing chamber we warmed the 
sample up, rotated it at room temperature, and cooled 
down again. The alignment uncertainty of the sample 
plane with the magnetic field was kept within 0.3°. The 
electron density as a function of gate voltage was deter- 
mined from Shubnikov-de Haas oscillations in perpendic- 
ular magnetic fields. We have verified that the gate volt- 
age dependence of the resistance in zero magnetic field is 
well reproducible in different runs with the accuracy of 
insignificant threshold shifts. In parallel magnetic fields, 
this dependence was used for determining the threshold 
voltage. 

In Fig.^ we show the dependence of the zero- field mo- 
bility, /X, on electron density, n., for two of our samples 
along with the data of Ref. jla| obtained on a Si/SiGe 
quantum well with higher mobility. We have verified that 
in the studied range of electron densities, the mobility 
is temperature- independent below 1 K. Despite the dif- 
ference between the mobilities is quite appreciable, the 
slopes of all dependences above 3 x 10^^ cm~^ are coin- 
cident and correspond to a power law: fi oc rig'^, as was 
observed previously [isl IT6j. Below 3 x 10^"'^ cm~^, the 
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FIG. 2; (a) The normalized magnetoresistance measured 
in different runs at ris = 2.98 x 10^^ cm~^ and at a tem- 
perature of 30 mK for B || 7 (solid and dotted lines) and 
B _L J (dashed line). The inset shows the magnetoresis- 
tance on another sample in one run at Us — 2.34 x 10^^ and 
2.54 X 10^^ cm~^ for B I. (b) The calculated magnetoresis- 
tance at Tis = 2 X 10^^ cm"'^ for different parts of the charged 
impurities at the interface. 



stronger decrease of the mobility with lowering electron 
density in our samples is likely to be a precursor of An- 
derson localization as caused by multiple scattering ■ 

In Fig.|2a), we show the magnetoresistance of the sam- 
ple in parallel fields. As the magnetic field is increased, 
the resistance increases and tends to saturate at the on- 
set of full spin polarization in this electron system. The 
magnetoresistance is practically independent of the rela- 
tive orientation of the field and the current, /. In other 
words, no anisotropy of the resistance with respect to the 
in-plane field is observed in our samples, which is similar 
to the case of Si MOSFETs. However, as compared to Si 
MOSFETs, the ratio R{Bc)/R{0) in the studied range of 
electron densities is much smaller, about 1.3 18]. This 
indicates directly that the remote-doping scattering pre- 
vails 01 . The form of the R{B) curve varies slightly in 
different runs. Even in the same run, it changes slightly 
with changing electron density (inset to Fig. [2J|. 

We now discuss the regime of remote-doping scat- 
tering. The screening properties of a 2D electron sys- 
tem are determined by two parameters: the screening 
wavevector, = 2gyme^ / eh? , and the Fermi wavevector, 
kp — {2'Kns/ g^y/'^ (where is the valley degeneracy, m 
is the band mass, and e is the dielectric constant). For 
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FIG. 3: Dependence of the integrand in Eq. (0 on trans- 
ferred wavevector for — 10^^ cm~^ (dashed line) and 
Us — 10^^ cm~^ (soUd hne) assuming G{q) = 0. 

our case the former equal to qs — 1.25 x 10^ cm^^ exceeds 
the latter which falls within the range 1.4 x 10^ cm"^ < 
2kF < 2.5 X 10^ cm^^ corresponding to the Wigner-Seitz 
radius 4.5 > Ts > 2.5. The third parameter that con- 
trols the electron scattering in Si/SiGe is half the inverse 
spacer width, l/2d, which is equal to 1.4 x 10^ cm~^ 
in our samples. As in the high-ris limit the transferred 
wavevector l/2d is small compared to 2kp, the electron 
backscattering is small. In the opposite limit the trans- 
ferred wavevector l/2d is approximately equal to 2kF 
and, therefore, the electron backscattering should occur. 
Nevertheless, we argue that this is still dominant small- 
angle scattering. Indeed, in the lowest order of multi- 
ple scattering theory lff| the inverse transport scattering 
time for the two- valley case can be written 

1 m f'^'"' q^dq 

r^^Jo kl^Akl-q^ [1 + ^(1-G(9))P' 

where G{q) is the LFC and in the case of remote 
doping is given by 

^'^-^ iV,(z)exp(-2gz)dz. (2) 
£ 9 Jd 

Here di — d is the width of the doped layer with doping 
density Ni. The dependence of the integrand, w{q), in 
Eq. ^ on transferred wavevector at high and low elec- 
tron densities is displayed in Fig. |31 ignoring the LFC. 
In both limits the main contribution to the scattering 
probability originates from wavevectors well below 2kF, 
in contrast to Si MOSFETs where the close vicinity of 
2kp contributes only. 

It is remarkable that in the regime of remote-doping 
scattering, the magnetoresistance at the onset of com- 
plete spin polarization is small Q . A qualitative account 
of this effect is given below. Assuming that the spin flip 
processes are absent, the calculated scattering rate 1/t 
for spin-up and spin-down electrons as a function of cor- 
responding electron density is shown in Fig. ^ With 
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FIG. 4: The calculated scattering rate for spin-up and spin- 
down electrons vs corresponding electron density at ris = 2 x 
10^^ cm~^ for different parts of the charged impurities at the 
interface. All dependences are normalized by the scattering 
rate for remote doping in zero magnetic field. 



increasing Uup (or ridown) the scattering angle decreases 
and, therefore, the scattering rate decreases except at the 
onset of complete spin polarization at which the increase 
of the scattering rate is related to the change of screening, 
similar to the effect discussed in Ref. ;4|- The resulting 
magnetoresistance is negative and only becomes positive 
near the ratio R{Bc) / R{0) being equal to 1.25, see 
Fig. |2b). It is important that if there is a small amount 
of charged impurities at the Si/SiGe interface, the neg- 
ative magnetoresistance is suppressed accompanied by 
somewhat larger R{Bc)/ R{0)- The calculated magne- 
toresistance with 0.5% of the charged impurities (about 
2 x 10^ cm^^) located at the interface describes the exper- 
iment reasonably well, see Fig. |21 So, density variations 
of the charged residual impurities at the Si/SiGe inter- 
face and, hence, of the local moments in the bandtail ^ 
naturally explain the observed changes in _R(_B) including 
those in different runs. 

If the local field corrections are still disregarded, the so- 
lution to Eq. is ^ (X nl^'^ as shown by the dashed 
line in Fig. ^ This dependence is noticeably weaker than 
the experimental one. The origin for the discrepancy is 
strong correlations in this electron system. At low elec- 
tron densities, worse screening is reflected by LFC so that 
the integrand in Eq. should be larger at wavevectors 
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about 2k F- The Hubbard form of the LFC 



G(<z) 



2gv y/(p 



(3) 



yields a power law oc n^ ®, the exponent being still 
smaller compared to the experimental finding. This is 
not very surprising because the Hubbard form includes 
only kinematic correlations that are caused by Pauli prin- 
ciple. Another contribution is given by dynamic correla- 
tions that are related to direct Coulomb interelectron in- 
teractions; these contribute to the LFC at large wavevec- 
tors j23|. To simulate the effect of dynamic correlations, 
we demand = 1 in Eq. (j^J. The so-obtained depen- 
dence fj,{ns) is indicated by the solid lines in Fig. ^ dis- 
regarding the numerical factor. Agreement between the 
calculation and the experiment shows that the LFC are 
approximately twice as large as the Hubbard form. 

Finally, we discuss the above assumption of the absence 
of spin flip processes. Available experimental data for Si 
MOSFETs allows one to conclude that such processes are 
in fact present but are not dominant. Indeed, the Hall re- 
sistance in a magnetic field with weak perpendicular and 
strong parallel components was found to be significantly 
lower than that expected for decoupled spin subbands 



|21|. At the same time, ratios R{Bc)/ R{0) « 4 were ob- 
served at low electron densities in the metallic regime ; 
these are close to the ratio R{Bc) / Rifi) — 4 expected for 
the case when no spin flip occurs, whereas twice as low 
a ratio is expected if the spin flip time is the shortest 
4]. This gives evidence that the spin fiip processes are 
not of importance and, thus, the above considerations 
are justified. 

In summary, we have established that in the 2D 
electron system in a modulation-doped Si/SiGe quan- 
tum well, the regime of remote-doping scattering oc- 
curs in which the parallel-field magnetoresistance is small 
and sensitive to uncontrollable density variations of the 
charged residual impurities at the Si/SiGe interface. 
Based on analysis of the mobility as a function of elec- 
tron density in this regime, we conclude that the local 
field corrections are approximately double the Hubbard 
form and, therefore, both kinematic and dynamic corre- 
lations make significant contribution. 

We gratefully acknowledge discussions with A. Gold, 
V. S. Khrapai, and S. V. Kravchenko. We would also like 
to thank D. Scheible, C. Meyer, and S. Manus for techni- 
cal help. This work was supported by A. von Humboldt 
Foundation via Forschungspreis, the Russian Foundation 
for Basic Research, and the Russian Ministry of Sciences. 



[1] V. T. Dolgopolov, G. V. Kravchenko, A. A. Shashkin, 
and S. V. Kravchenko, JETP Lett. 55, 733 (1992); 
D. Simonian, S. V. Kravchenko, M. P. Sarachik, and 
V. M. Pudalov, Phys. Rev. Lett. 79, 2304 (1997); S. V. 
Kravchenko, D. Simonian, M. P. Sarachik, A. D. Kent, 
and V. M. Pudalov, Phys. Rev. B 58, 3553 (1998); K. M. 
Mertes, D. Simonian, M. P. Sarachik, S. V. Kravchenko, 
and T. M. Klapwijk, Phys. Rev. B 60, R5093 (1999); 
V. M. Pudalov, G. Brunthaler, A. Prinz, and G. Bauer, 
JETP Lett. 65, 887 (1997); A. A. Shashkin, S. V. 
Kravchenko, and T. M. Klapwijk, Phys. Rev. Lett. 87, 
266402 (2001). 

[2] T. Okamoto, K. Hosoya, S. Kawaji, and A. Yagi, Phys. 
Rev. Lett. 82, 3875 (1999); S. A. Vitkalov, H. Zheng, 
K. M. Mertes, M. P. Sarachik, and T. M. Klapwijk, Phys. 
Rev. Lett. 85, 2164 (2000). 

[3] At electron densities below approximately l.Sric (where 
ric is the density for the zero-field metal-insulator transi- 
tion), the application of a strong parallel magnetic field 
makes the system insulating, in which case the ratio 
R{Bc)/ R{0) may become arbitrarily high l]. 

[4] V. T. Dolgopolov and A. Gold, JETP Lett. 71, 27 (2000). 

[5] A. A. Shashkin, S. V. Kravchenko, V. T. Dolgopolov, 
and T. M. Klapwijk, Phys. Rev. Lett. 87, 086801 (2001); 
S. V. Kravchenko, A. A. Shashkin, and V. T. Dolgopolov, 
Phys. Rev. Lett. 89, 219701 (2002). 

[6] V. M. Pudalov, G. Brunthaler, A. Prinz, and G. Bauer, 
Phys. Rev. Lett. 88, 076401 (2002); V. T. Dolgopolov 
and A. Gold, Phys. Rev. Lett. 89, 129701 (2002); A. 
Gold and V. T. Dolgopolov, J. Phys. Condens. Matter 
14, 7091 (2002). 



[7] J. M. Broto, M. Goiran, H. Rakoto, A. Gold, and V. T. 

Dolgopolov, Phys. Rev. B 67, 161304(R) (2003). 
[8] J. Zhu, H. L. Stormer, L. N. Pfeiffer, K. W. Baldwin, and 

K. W. West, Phys. Rev. Lett. 90, 056805 (2003). 
[9] E. Tutuc, S. Melinte, E. P. De Poortere, M. Shayegan, 

and R. Winkler, Phys. Rev. B 67, 241309 (2003). 
[10] M. Y. Simmons, A. R. Hamilton, M. Pepper, E. H. Lin- 
field, P. D. Rose, D. A. Ritchie, A. K. Savchenko, and 

T. G. Griffiths, Phys. Rev. Lett. 80, 1292 (1998); J. 

Yoon, C. C. Li, D. Shahar, D. C. Tsui, and M. Shayegan, 

Phys. Rev. Lett. 84, 4421 (2000); S. J. Papadakis, E. P. 

De Poortere, M. Shayegan, and R. Winkler, Phys. Rev. 

Lett. 84, 5592 (2000). 
[11] T. Ando, A. B. Fowler, and F. Stern, Rev. Mod. Phys. 

54, 473 (1982). 
[12] A. Gold, Appl. Phys. Lett. 54, 2101 (1989). 
[13] Y. Y. Proskuryakov, A. K. Savchenko, S. S. Safonov, L. 

Li, M. Pepper, M. Y. Simmons, D. A. Ritchie, E. H. 

Linfield, and Z. D. Kvon, J. Phys. A: Math. Gen. 36, 

9249 (2003). 
[14] A. Gold, Physica E 17, 305 (2003). 

[15] Z. Wilamowski, N. Sanderfeld, W. Janish, D. Toebben, 
and F. Schaeffler, Phys. Rev. Lett. 87, 026401 (2001). 

[16] An attempt to account for the strong dependence fj,{ns) 
was made in Ref. [isl based on the idea of worse screen- 
ing of potential fluctuations with lowering due to 
a decrease of the average thermodynamic density of 
states near the metal-insulator transition. However, as 
the range of potential fluctuations found in that paper ex- 
ceeds the electron mean free path, the scattering calcula- 
tions using the average thermodynamic density of states 



5 



in such an inhomogeneous electron system are obviously 
incorrect. 

[17] A. Gold, Phys. Rev. B 44, 8818 (1991). 
[18] A similar dependence R{B) was observed by T. Okamoto, 
M. Ooya, K. Hosoya, and S. Kawaji, _cond-mat/0307521 



[19] A. Gold and W. Gotze, Phys. Rev. B 33, 2495 (1986). 
[20] B. Davoudi, M. Polini, G. F. Guhani, and M. P. Tosi, 

Phys. Rev. B 64, 153101 (2001). 
[21] S. A. Vitkalov, Phys. Rev. B 64, 195336 (2001). 



